ANALYTICAL APPROXIMATION OF BLASIUS' SIMILARITY 
SOLUTION WITH RIGOROUS ERROR BOUNDS 
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Abstract. We use a recently developed method [2], [3] to find accurate an- 
alytic approximations with rigorous error bounds for the classic similarity so- 
lution of Blasius of the boundary layer equation in fluid mechanics, the two 
point boundary value problem /"' -|- //" = with /(O) = /'(O) = and 
'. lim 

a;_>oo /'(^) — 1- The approximation is given in terms of a polynomial in 
[0, ^] and in terms of the error function in oo). The two representations for 
' the solution in different domains match at a; = ^ determining all free parame- 

CN I ters in the problem, in particular /"(O) = 0.469600 ± 0.000022 at the wall The 

■ method can in principle provide approximations to any desired accuracy for 

this or wide classes of linear or nonlinear differential equations with initial or 
boundary value conditions. The analysis relies on controlling the errors in the 
approximation through contraction mapping arguments, using energy bounds 

■ for the Green's function of the linearized problem. 
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1. Introduction and main results 



Finding exact expressions for solutions to nonlinear problems is an important 
area of research. Closed form solutions however exist only for a small sub-class 
of problems (essentially for integrable models). On the other- hand, if a problem 
involves some small parameter e (or a large parameter) and the limiting problem 
is exactly solvable, then there exist quite general asymptotic methods to obtain 
convenient expansions for the perturbed problem. 

Indeed, consider for instance the question of finding the solution to M[u, e] ~ 0, 
where A/" is a (possibly nonlinear) differential operator in some space of functions 
I satisfying boundary/initial conditions and that uq is the solution at e = 0. Exis- 

■ fence and uniqueness of a solution u as well bounds on the error E = u — Uq may 

CO ' be found as follows. We write 

LE = -5-Ni{E) 

; where L = |^|„=„o, ^ = A^^] and Ni{E) ^ N {u^ + E) -LE = 0{e^), inverting 

?— ( ' i in a suitable way, subject to the initial/boundary conditions, and using the 

contractive mapping theorem in an adapted norm to control the small nonlinearity 
Ni{E). 

A relatively general strategy has been recently been employed |2]-[3] in prob- 
lems without explicit small parameters. The approach uses exponential asymptotic 
methods and classical orthogonal polynomial techniques to find a function uq which 
is a very accurate global approximation of the sought solution u, in the sense that 
M{uq) is very small in a suitable norm and the boundary conditions are satisfied 
up to small errors. Once this is accomplished, a perturbative approach similar with 
the one above applies with the role of e played by the norm of 5 and one obtains 
an actual solution u by controlling the equation satisfied hy E = u — uq. 

In the present paper, we apply this strategy to the well-known Blasius similarity 
equation arising in boundary layer fluid-fiow past a fiat plate. We improve the 
methods in [2]- [2] by replacing the laborious detailed analysis of the Green's function 
with softer energy methods. 
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The Blasius similarity solution satisfies 

(1) /'"(x) + /(x)/" (x) = for X e (0, oo) 
with boundary conditions 

(2) /(0) = 0,/'(0) = 0, hin /'(.t) = 1 

2:— ^+00 

This equation and various modifications have garnered much attention since Blasius 



[4] derived i i^-^'l Existence and uniqueness were first proved by Weyl in 0- Issues 
of existence and uniqueness for this and related equations have been considered 
as weU by many authors (see for eg. [7], [8], the latter being a review paper). 
Hodograph transformations |S] allow a convergent power series representation in 
the entire domain, but the convergence is slow at the edge of the domain and the 
representation is not quite convenient in finding an approximation to / directly. 
Empirically, there has been quite a bit of interest in obtaining simple expressions for 
Blasius and related similarity solutions. Liao [9] for instance introduced a formal 
method for an emprically accurate approximation; the theoretical basis for this 
procedure and its limitations remain however unclear. We are unaware of any 
rigorous error control for this or any other efficient approximation in terms of simple 
functions. 

In [5] , using a transformation introduced by Topfer , it is also proved that / 
in dH), ([2) can be expressed as 

(3) f{x)^a-'/'F[a-'/'x), 

where F satisfies the initial value problem 

(4) F"'ix) + F{x)F"{x) =0 for X e (0, 00) 
with initial conditions 

(5) F(0) ^ , F'(0) = , F"{Qi) = 1, 

In 0, lim^^ooi^'(x) = a e R+ (cf. [3 [U]). Conversely, if /(x) satisfies ^- 
^ with /"(O) = (5 > Q (physically, this corresponds to positive waU stress), it 
may be checked that F[x) = /3~^/^/(/?~^/^x) satisfies (U)-®. Because of this 
equivalence, it is more convenient to find an approximate analytical representation 
for the solution F of (H])-® and determine the value of a. The solution / to the 
original problem is obtained through the transformation ([3]); the stress at the wall 
is given by 



(6) /"(0)=a 



-3/2 



1.1. Approximate representation of the solution. Let 



equation in the original Blasius' paper has a eoefHcient i for //'; however the change 
of variable x — > x/\/2, f — //\/2 transforms JlJ into Blasius' original equation. Thus, /"(O) = 
0.469600 ± 0.000022 transforms to /"(O) = 0.3320574 ± 0.000016 in the original variables. 
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where [po, ...,^12] are given by 

510 18523 42998 113448 65173 390101 



(8) 



2326169 



10445149' 5934' 441819 ' 81151 ' 22093' 6016 ' 9858 ' 

4134879 1928001 20880183 1572554 1546782 1315241 1 



7249 



1960 



19117 



2161 



5833 



32239 J 



Define 

(9) t{x) = ^(x + &/a)^ loit) = l-\^e*erfc(\/i) , Jo(<) = l-\/2^e^*erfc(y2l) 
where erfc denotes the complementary error function and let 

(10) go(i) = 2cVie-*h + c'e.-'* (2 Jo - h - Q , 

The theorem below provides an accurate representation of F of ([4]), ([5]). 
Theorem 1. Let Fq be defined by 



(11) 



Foix) 



4 + x^P (fx) for X e [0, f ] 
2iX5) 90(^(2;)) for X > 



Then, there is a unique triple {a,b,c) close to {ao,bo,co) = (fffgi ~f^7 ^^^) ^'^ 

the sense that (a, b,c) € S where 

(12) 



(a,6,c) e 



(a-ao)2 + i(6-6o)2 + J(c-co)2 < po 



5 X 10' 



with the property that Fq is a representation of the actual solution F to the initial 
value problem ^-(^ within small errors. More precisely, 

(13) F{x) = Foix) + E{x) , 

where the error term E satisfies 

5, 



(14) \\E"\U 
and for x > ^ 



< 3.5 X 10-'' , llS'lloo < 4.5 X 10- 



l^lloo < 4 X lO"*^ on [0, 



(15) 



E 



< 1.69 X 10"^^ 



dx 



-E 



< 9.20 X IQ-H-^^^e-^* 
d^ 



dx^ 



:E 



< 5.02 X lO^-'r'e" 



Remark 1. Certainly, F is smooth since it is is an actual solution of Q,©, which 
exists on [0,oo) and is unique, see [5]. However, the particular choice (a, 6,c) £ S 
in Theorem [T] needed in order ioi F = Fq + E to solve (HI)-® docs not ensure 
continuity of the approximate solution Fg at x = |. Nonetheless, if Fq,Fq,Fq 
are needed to be continuous, this can be achieved by a slightly different choice of 
{a,b,c) £ S (see Remark [T^ . namely 

(16) (a, b, c) = (1.6551904561499..., -1.565439826457.., 0.233728727537...) . 



Note also that ([T5|) implies not only small absolute errors (that, in the far field 
hold even for the approximation of F" by zero) but also very small relative errors 
for a; > 5/2. 
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Definition 1. Let ai = gq - po, = ao + po, bi = bo — 2po, 6r = &o + 2po, 
Q = Co — 2po and c,- = cq + 2po- We see that (a, 6, c) e S, implies that a e [a/, Or], 
6 e [6;,&r], c e (Q,Cr). We define = f (I + I) i^o^c that a; e [|,oo) 









2 






a; 












2 


2 




' 2 







corresponds to t e [t„i, cx)) and t„i G 
(1.998859- •• , 1.999438- ••)■ 

Remark 2. The error bounds proved for E in Theorem [1] are hkely a 10 fold over- 
estimate. Comparison with the numerically calculated F suggests that \F — Fo\, 
\F' - F(i\ and |F" - F^'\ in [O, |] are at most 2 x 10-^2 x lO"'^ and 5 x 10"^ 
respectively. Using the nonrigorous bounds on E and its derivatives reduces the po 
in the definition of S from 5 x 10~^ to 1.4 x 10~^. It is thus likely that (a, 6, c) sa 
(ao,&OjCo) with five (rather than the proven four) digits accuracy. 

The proof of Theorem [1] rests on the following three propositions, proved later 
in the paper. 

Proposition 2. The error term E{x) = F{x) — Fq{x) verifies the equation 

(17) C[E] := E" - FqE" - F^'E = -F^" - FqE^' - EE" 

(18) E{0) = = E'{0) = E"{0) 
and satisfies the bounds 1^14^ on I = [O, |] 

Proposition 3. For given {a,b,c) with a > 0, |c| < j. in the domain x > — ^ + 

\J^-, for T > 1.99, which corresponds to t > T > 1.99, there exists unique solution 
to ^ in the form 

€ + b 



(19) 



Fix) 



2t{x) 



qitix)) 



that satisfies the condition lim4_j.oo — > 0. Furthermore, 
(20) q{t)=qo{t)+£it) 
where £ is small and satisfies the following error bounds: 



(21) 
(22) 
(23) 



-3t 



< 1.6667 X 10~ 



9^3/2 



£'{t)--£{t) 



< 1.6667 X 10" 



Vi£"{t)-l=£'{t) + ^£{t) 



< 1.6667 X lO-^'t-^e-^* 



Proposition 4. There exists a unique triple (a, b,c) €z S so that the functions in 
the previous two propositions: Fq{x) + E{x) for x < ^ and ax + b + y^^^^y 9(^(2;)) 
for x > I and their first two derivatives agree at x = ^. 

The proof of Theorem [T] follows from Propositions in the following manner: 
Proposition [5] implies F{x) = Fo{x) + E{x) satisfies (H])-® for x £ I; we note that 
Fo{0) = = F^{0) and F^'(O) = 1. 
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Proposition [3] implies F{x) ^ ax + b+ y 2t|j) llo{tix)) +£{t{x))] satisfies ^ in 

a range of x that includes [|, oo) when (a, b, c) € S. Further, Proposition |4] ensures 
that this is the same solution of the ODE (|4]) as the one in Proposition^ Identifying 

Fo{x) and E{x) in Theorem [T] in this range of x with ax + b + ^ 2t1x) 1oi^ix)) ^^'i 
\J 2t{x) ^^^^^))^ respectively, and relating x-derivatives to t derivatives, the error 
bounds for E, E' and E" follow from the ones given for £ in Proposition 2] for 
(a, 6,c) g S. The proofs of Propositions [2H] are given in Sections [2E] respectively. 

2. Solution in the interval / = [0, |] and proof of Proposition [2] 

The ansatz for Fq in the compact set / is obtained simply by projecting an em- 
pirically obtained high accuracy approximate solution on the subspace spanned by 
the first few Chebyshev polynomials. More precisely, to avoid estimating deriva- 
tives of an approximation, which are not well-controlled, we project instead the 
approximate third derivative F'" = —F"F on the interval^ / = [0, |]. The rig- 
orous control of the errors of the integrals of F'" is a much simpler task, and this 
is the way we prove Proposition [21 For a given polynomial degree, a Chebyshev 
polynomial approximation of a function is known to be, typically, the most accurate 
in polynomial representation, in the sense of L°° . 

We seek to control the error term E in (jl3|) by first estimating the remainder 

(24) R{x)=F;,"{x)+Fo{x)F^'{x), 

which will be shown to be small (< 0.673 x 10"^). Then, we invert the principal part 
of the linear part of the equation for the error term E by using initial conditions 
to obtain a nonlinear integral equation. The smallness of R and careful bounds on 
the resolvent allow for using a contraction mapping argument to obtain the sharp 
estimates for E and its derivatives stated in Proposition [2] 

2.1. Estimating size of remainder R{x) for x £ I. Since P is a polynomial 
of degree twelve, R{x) is a polynomial of degree 30. We estimate i? in / in the 
following manner. We break up the interval into subintervals {[ ]} with 

Xq = and X14 = |, while {.Tj}j|^-|^ is given by 

{0.0625, 0.125, 0.25, 0.375, 0.50, 0.75, 1.0, Xc, 1.5, 1.75, 2.0, 2.25, 2.40} 

where Xc = 1. 322040 [£] The intervals were chosen based on how rapidly the poly- 
nomial R(x) varies locally. 

Wc re-expand R{x) as polynomial in the scaled variable r, where x = ^ {xj + a;j_i)-|- 
^ {xj — Xj-i) T. and write 

30 
fc=4 



(■^^Chcbyshev polynomial approximation is on [—1, 1] interval, but a linear scaling and shift of 
independent variables can accommodate any finite interval 

will be found later, it is convenient to choose one of the subdivision points Xc to be 
approximately, to the number of digits quoted, the value of x where Fq{x) — 2Fq(x) + 1 changes 
sign. 
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and determine the maximum Mj and minimum mj of the third degree polynomial 

P^-'\t) for r e [—1,1] (using simple calculus). We estimate the error on the 
remaining coefficients: 

30 

fe=4 

It follows that in the j-subintcrval we have 

mj - e'-jI^ < R{x) < Mj + E'jp 

The maximum and minimum over any union of subintervals is found simply taking 
min and max of rrij — E^jp and Alj + E^^"^ over the the indices j for subintervals 
involved. This elementary though tedious calculatiorF^ yields 

(25) - 3.22 X 10"'^ < R{x) < 2.505 x 10"'^ for x £ [Q,Xc] 

4.6 X 10"** < R{x) < 4.06 x 10"^ for x € [xc,2.0] 

2.78 X lO"'^ < R{x) < 6.73 x lO"'^ for x £ [2.0,2.5] 

We note that the remainder is at most 6.73 x 10"^ in absolute value in the interval /. 
In a similar manner of finding maximum and minimum, bounds for the polynomials 
Fo{x), Fq and Fq{x) may also be found. For x £ [O, i] , 

(26) - 5 x 10^1° < Fo{x) < 0.008 , -8 x lO"^^ < ^^(^,) < g.ia , 

0.99 < i^o'(a;) < l + 2x 10"^ 

while for x £ [|, |] , 

(27) 0.03 < Fo{x) < 2.59 ,0.12 < F^{x) < 1.7 ,0.09 < F^^'ix) < 1 

2.2. Properties of some functions used in the sequel. Based on the calcu- 
lations above, one can also conclude that Fq — 2Fq + 1 and 2Fq — 2Fq have only 
one zero in the interval / in the following manner. Note that derivatives of these 
functions are -FqE^' - 2F^ + R < and -2FoF^' - 2F^ + 2i? < respectively in 
|, |] , where R has been estimated in the previous subsection. From the bounds 
in the interval [O, |] , it is clear that Fq' — 2Fo + 1 and 2Fq' — 2Fo are positive in 
[0, i] . Thus, we conclude that F^' - 2Fo + 1 and 2F^' - 2Fo can have at most 
one zero in /. The values of F^' - 2Fo + 1 at Xc = 1.322040 and 1.322041 have 
opposite signs, implying that there is a unique zero in / in between two numbers 
(recalling that the derivative is negative). Similarly, we conclude there is a unique 
zero of 2F,5' - 2Fo between 1.2314283 and 1.2314284. Similar arguments show that 
F^'{x) - 2Foix) only has one zero in J at a; = 0.9399325 ■ ■ ■ < Xc- 



(4)The maximum and minimum found through analysis described here is found to be consis- 
tent with a numerical plot of the graph of R{x), as must be the case. The calculations can be 
conveniently done with a computer algebra program, as they only involve operations with rational 
numbers. 
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2.3. Green's function estimate for x G [a::;,a;r]. Consider now the problem of 
solving the linear generally inhomogeneous equation 

(28) £[0] + FoixWix) + F[,\x)cj){x) = r{x) 

over a typical subintcrval [x/jX^] C /, with initial conditions 4>{xi), (t>'{xi) and 
(f)" {xi) known. The solution of this inhomogeneous equation is given by the standard 
variation of parameter formula: 

(29) 0(a;)=^ 0(^-1) (xO<I>,(x) + ^$^(x) [\,{t)r{t)dt 

where form a fundamental set of solutions to Ccf) = and {^j"(2;)}^^j^ are 

elements of the inverse of the fundamental matrix constructed from the $j and 
their derivatives. The precise expressions are unimportant in the ensuing: we only 



(5) 



need their smoothness in x. It also follows from the properties of and ^ j 
that 

(30) cp"{x) = ^0(^-1) (xO$;'(x) + ^<l>;'(x) [\,{t)r{t)dt , 

3 = 1 j=l -^^^ 

It is useful to write (PO]) in the following abstract form 

3 

(31) r{x) = <^^'~^^ {xiW3{x) + Q W (a;) 

where from general properties of fundamental matrix and its inverse for the linear 
ODEs with smooth (in this case polynomial) coefficients is a bounded linear 
operator on C([x/, x,]); denote its norm by M, 

(32) M=\\g\\ 

Then, on the interval we have, 

oo 

(33) ||(/."|loo<^A'/jU'^"'^(a;0 +^^lklloo; Mj= sup $;'(x) 



We will now estimate Mj for j = 1..3 and M indirectly, using "energy" bounds. 
Because of linearity of the problem, for the purposes of determining these bounds, 
it is enough to separately consider the cases (i)-(iii), when r = 0, (\)^^^^^(xi) =■ for 
1 < 7^ j ^ 3i (\>^^~^^ = 1 respectively, and, finally, (iv) when ^'^^~^^{xi) — 
for fc 1, ..3 and r{t) ^ 0. 

For all cases (i)-(iv), we return to the ODE 

(34) 0"' + Fo0" + F'^4> = r 

Multiplying by 2(/)", integrating from xi to x and using initial conditions, it follows 

that 

(35) 

(^"(x))' = (0"(xO)' - r {2F,{y) {4l\y)f + 2F'^{y)4l\y)4>[y) - 20"(y)r(y)} dy. 



(5)ln particular, Y,]=\ = 0, Ej=i = 
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We note further that, given 0(a:;) and (j)'{xi), 4'{x) is determined from (/)"(a;) and 
the relation 



(36) ^{x):=<i}{x)-ct){xi)^{x-xiW{xi)= I {x - y)r{y)dy 

Using ([361) in ([351), it follows that 



(37) {r{x)r = {r{xi)r 



2F^'(y) [cj^ixi) + {y - xi)^'{xi)] cj," {y)dy 



{2Fo(y) {ct>"{y)f + 2F^'(y)0"(y)0(y) - 24>" {y)r{y)] dy, 



2.4. Case (i): Determination of Mi. In this case, we set (l){xi) = 1, (t)'{xi) 
4>"{xi) = = r{x) in ([57)1 to obtain (using Fq > and Cauchy-Schwartz): 



(38) {c^"{x)Y< Fi;{y)dy+ F^^ {y)4>\y)dy 



It is convenient to define, for xi < t < Xr, 



(39) 



D{t) = sup 

xe[xi,T] 



iriy)f{2F;,'iy)-2Foiy)}dy 



(x) 



From (pG]) and the definition of D we get 

(40) <VdW). ~m < 



VD{xj, (0"(x))' < -D(x); in [xi,Xr], 



Using (gOl) and ([381) ^e see that 



(41) i^"ix)r< K(y)rfy + F;,'iy) 



D{y)dy 



Define now 



(42) gi(x) = F^'(x) 2 + 



(a; - xiY 



and Qi{x) 



{r{y)f{2F;,'{y)-2Fo{y)}dy 



- 2Fo{x) if 2F^'(x) - 2Fo(a;) > 

= (^~/')V -'(.T) if 2F:,'ix) ~ 2F,{x) < 



Equations (j40p and (j4ip imply the following inequality for for < < x < r < 



(43) {cP"{x)r< F(,'{y)dy+ D{y)Q^{y)dy 

Jxi Jxi 

Since the right side is independent of x, it follows that 

(44) D{t)< [ F;^{y)dy+ [ D{y)Q,{y)dy 



i{v)dy 



From Gronwall's lemma, it follows that 
(45) 



D{t) <([ Fi;{y)dy] cxp / Q 

\Jxi J VJ xi 
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Evaluating (|45p at t ^ Xr immediately implies 



(46) 



Qi{y)dy 



2.5. Case (ii): Determination of M2. In this case, we set 4>'{xi) = 1, (j){xi) 
= 4>"{xi) = = 7'(a:;) in (|37|) to obtain by Cauchy-Schwartz 



(47) i^"{x)r< Fi;{y){y-xiYdy\+ Fll{y)4>'{y)dy 



{4>"{y)r{2Fi;{y)-2Fo{y)}dy 
Again introducing D{x) and Qi{x) as in case (i), we obtain the inequality 
(48) D{t)<[ {y-xifF(,'{y)dy+ f D{y)Qi{y)dy 



Qi{y)dy 



and therefore, it follows from Gronwall's Lemma that 



(49) 



M2 < 



{y-xi)^F[;{y)dy 



1/2 



exp 



2.6. Case (iii): Determination of il/3. In this case, we set 4>"{xi) = 1, (j){xi) 
^ (f)' (xi) ^ = r{x) in ([321) to obtain 



'{y)r{2F^;{y)-2Fo{y)}dy 



(50) (0"(x))^<l+ F^'(y)0^(y)dy + 



Once again introducing D as in case (i) and defining 

(x-xi)"^ 



(51) Q2ix) 

([5n|) imply 
(52) 



and Q2{x) 



F^'(x) - 2Fo{x) , if F^'(x) - 2Fo{x) > 
(x-xiY 



-F^'ix) , ii F^'{x)~2Fo{x)<0 , 



D{t) < 1+ / D{y)Q2{y)dy 



Gronwall's Lemma and definition of D implies 

^1 



(53) 



Ms < exp 



'■{y)dy 



2.7. Case (iv): Determination of M = \\g\\. With (j){xi) = = (j)'{xi) = 
(j>"{xi) = 0, (|57)) implies by Cauchy-Schwartz 



(54) {cjy"{x)Y< r\y)dy+ F^,' iy)4>\y)dy 



'{y)y{F;,\y)-2Fo{y) + l}dy 
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We define D as in (i) and also 

(55) Q{x) = Fi;{x) - 2Fo{x) + 1 + t^-p^F;,'ix) if F^,' - 2i^o + 1 > 



The inequality 
(56) 



Q{x) 



r^{y)dy 



{x - xiY- 



F^{x) if i^^' - 2Fo + 1< 



J xi 

)lying from < 
\/D{xr) < ||r||oo (xr ~ xiY^^ exp 



D{y)Q{y)dy 



follows from (j54p . implying from Gronwall's Lemma 
(57) ' ^ 

implying 
(58) 



Q(y)dy 



M < {xr — xi)^^'^ exp 



Q{v)dy 



2.8. Existence of _F; error bounds for x G [a;;,a;r] C /. Consider the decompo- 
sition 

(59) F{x) ^ Fq{x) + E{x) 

We seek to find error estimates for Fix) and its first two derivatives for x £ F 
For this purpose we break up / into a number of subintervals. Note that for the 
first subinterval xi — 0, where F{xi) = = F'{xi) = F"{xi). Consider a typical 
subinterval / = where the bounds on F{xi), F'{xi) and E" ixi) on earlier 

subintervals have been already obtained. The equation for F{x) on [x/,a;r] is 

(60) C[E] := E"'{x) + Fq{x)E" {x) + F;^'{x)E{x) ^ ~E{x)E"{x) - R{x) 

Inverting C as described in previous subsection leads to the following integral equa- 
tion: 

3 

(61) E"{x) = ^i?(^-i)(xz)$;'(x) - g [R] (x) + g [EE"] (x) =:Af[E"] (x) 

i=i 

and where E is given by 



^ (Xr Xl) .^^11 



(62) E{x) - E{xi) - {x - xi)E' {xi) E{x) ^ f {x - t)E" {t)dt 
Note that ([62]) implies 

(63) ll^ll 



We prove the following Lemma that, once some bounds are satisfied, ensures the ex- 
istence, uniqueness and smoothness of the solution E of (|60p and provides estimates 
of -B, E' and E" . 



Lemma 5. Assume that for some e > we have 

(64) M {\^E{xi)\+ {xr ~ xi)\e' {xi)"^ (1 + e) + i(x, - x,)'MBo(1 + e)' 



< 



(65) 



M 



{\^F{xi)\ + {xr-xi)\F'{xi)'^+{xr-xiYMBQ{l + e) 



< 1 
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where 

3 

(66) Bo = M\\R\\^ + J2m,\e^^\xi) . 

Then, there exists a unique solution E" of (pT]) in a baU of radius Bo{l + e) in the 
sup norm in C{[xi,Xr])- 

Under these assumptions, E is in C'^([xi,Xr]) and satisfies Ii60\) with initial con- 
ditions E^^\xi), j = 0, 1, 2 and 

(67) ||S"'||oo < ||i^o||oo(l + e)Ba + \\F^'\\oo (\Eixi)\ + (x, - a;,) | ^' (x, ) | ) 

+ ^!I-F^o'lloo(a;r - xi)^Bo{l + e) + ^{xr - xi)^B^o{l + e)^ + \\R\\o. 

Proof. Since Q is the Green's function of a smooth hnear ODE, Q maps C{[xi, Xr]) 
into itself; the same, clearly, holds for M. From the definitions of M in ((32)) . and 
of Mj, j = 1, 2, 3 in ^ (whose bounds will be obtained using ge]), (gg), ^ and 
it follows that 



\\Af[E"]\\^<M\\R\\^ + Y,M,\E<^^-'\xi] 



+ M {\E{xi)\ + \Xr - xi\\E'{xi)\) ||i?"||oo + ^M||i?"|lL 



and 

(69) \\Af[E"]-Af 



E" 



< M {\E{xi)\ + la;, - xi\\E' {xi)\) \\E" - E"\ 

{Xr - XiY 



-M HE"' 



\E" 



E" -E"\ 



Using (|63p , ([64]) and (|66)) in ([68)) and (|69)) we sec that M is contractive in a ball of 
radius (1 + e)i3o in ^^([^c;, a:,]), implying existence and uniqueness of a solution to 
dSH). Clearly, is equivalent to ^] from ^ it follows that is also con- 
tinuous. Now, E'" is easily estimated from (|60p in terms of lower order derivatives, 
and the result follows. | 

2.9. Determining E using Lemma [5l In this section, we break the interval 
[0,5/2] in a suitable way and show that Lemma [5] applies in all subintervals. The 
choice of subintervals is Ii = [0, x^ Ti = [xc, 2], I3 = [2, |] , where Xc = 1.322040 
is, within number of digits quoted, the zero of F^f — 2Fq + 1 obtained in 



2.9.1. Error estimates on X\. On Ii, using (|11|). it is easily checked that 

(70) M < 3.03 

while by (f25|) we have ||i?||oo.ii < 3.22 x 10~^ Since the initial conditions on this 
interval are E{xi) = E'{xi) = E"{xi) = 0, the Mj do not contribute to ([Ml), ([69]) . 
and (|66p implies 

(71) Bo < 0.9757 x 10"^ 

The conditions (|64p and ([55)) are satisfied for e = 3 x 10^^ so that Lemma [5] implies 



(72) 



\E"\\ooM < 0.976 X 10 



-6 



12 



O. COSTIN^ AND S. TANVEER^ 



On integration it follows that 

(73) P'l|oo,ii < 1.29 X 10-*^ , \\E\\^,i, < 0.853 x 10"*^ 

2.9.2. Error estimates on l2- On I2 = [a^'c, 2], using (jll[) . it is easily checked that 

(74) Ml < 0.572 , A/2 = 0.199 , Ah < 1.01, M < 0.825. 

Since at xi, we can apply ([72|) and ([73|) to bound E{xi),E'{xi),E"{xi), using 
||-R||oo,i2 < 4.06 x 10-^, JMl) implies 

(75) Ba < 2.0653 x 10"^ 

Lemma [5] applies if e = 2 x 10^®. Therefore, the solution E exists and is unique on 
I2 and 

(76) < 2.07X 10-*^ 

By integration and using the bounds on E(xi), E'{xi) obtained in the previous 
interval, see ([TS]) . we get 



(77) 



\E'\\oo,io < 2.7 X 10-*^ , ||^||oo,i2 < 2.21 X 10- 



2.9.3. Error estimates on I3 = [2, 2.5]. On I3 = [2, 2.5], using ^ we get 

(78) Ml < 0.3 , A/2 < 0.0744 , M3 < 1.01 , A/ < 0.708 

while |j/?||oo,[2,2.5] < 0.673 x 10^® by (j^S]) . Proceeding as in the previous interval 
we get 

(79) Bo < 3.431 X 10^^ 
Here Lemma [5] applies with e = 3 X 10"*^ and thus 



(80) 



|S"||ooT, < 3.44 X 10- 



Proceeding as in the previous intervals, we get 



(81) 



|^'||ooj3 < 4.42 X 10-6 , ||S|U,i3 < 3.99 x 10" 



2.10. End of proof of Proposition [2l In the previous subsection we have shown 
that Lemma [S] applies on each of the intervals Xj entailing the existence and unique- 
ness of a solution E of the initial value problem (|17p - (fT8|) over the interval /. The 
same calculations show the L°° norms of E^^") satisfy the bounds in Theorem [TJ 

3. Solution in t > T > 1.99 for \c\ < a> and proof of Proposition [3] 
Wc decompose 

(82) F{x) =ax + b + g{x) 
Then, it is clear from (|4]) that g satisfies 

(83) 5"' + {ax + b + g) g" = 



Lemma 6. Assume a > 0. Then any solution g to 
X — > +00 has the following asymptotic behavior 



for which 17 — > as 



(84) 
(85) 



g (x) ~ Cexp 
C 



a 

-- \ x + - 



a 



{ax + b) 



exp 
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(86) 



c 



{ax + by 



■ exp 



- a: + - 



Proof. Eq. ((83)) implies 

(87) g"{x) = Cexp 













J Xq 



Since g{x) — o(l), f^^g{t)dt = o(x) as a; — >■ oo; a > and (f87|) imply that for 
sufficiently large x we have 















< 


C 


exp 















- ex 



Integrating g" and using ([88]) . we get for large x, 
for some Ci. Similarly, for large x we get 



■ exp 



\ x + - \ + ex 



(90) 



g{x) - Cix - C2 



< 



C 



{ax + b — e)'- 



■ exp 



- \ x + - 



ex 



Since g ~ o(l), we must have Ci =0,(72 = 0, giving rise to an exponentially 
decaying a priori bounds on g in (|90|) (with Ci = = C2). We can then set 
ccq = 00 in ([57)) . and (|M)) follows, and by integration, ([55]) and hold. | 

Lemma 7. For given a > 0, b and C , for xq > sufficiently large, there exists 
unique solution to (|83p in [xo,oo) with asymptotic behavior ^8^-(8^. 

Proof. Take xo > sufficiently large so that 



(91) 



-a/2{xo+b/a}- 



< 



{axo+b)3 A\C\ 



We define H{x) = exp [f (a;+ ^)^] g"{x). Then, from Lemma [SI we see that the 
appropriate space for H is C([a:o,oo)) with the sup norm. Lcmma[n]also implies 



(92) g{x) 
Eq. (j92p immediately implies 



exp 



00 00 



2 a 



H{s)dsdy 



-f (x+fc/a)2 

(93) \g{x)\ < ||g||oo 

{ax + by 

From (|83l) we obtain the following integral equation: 

(94) H{x) = C- / H{s)g{s)dsdy =:Af[H]{x) , 

where g is determined from using We analyze in C([a;o,oo)). We 

consider the space of continuous functions in [a;o,oo). equipped with ||.||oo norm. 
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Clearly M maps a ball of radius 2\C\ in this space back to itself since ([93]), ([94]) 
and ([9T|) implies 

(95) mmWoo < \C\ + ^-—^\\H\\l < 2\C\ 

Thus J\f maps the ball Bq of radius 2\C\ into itself and is contractive there since 

-a/2(xo+b/af 

(96) WNIH,] -Af[H,]\U < -Tv^ (Ili^ilU + ||i?2||oo) \\H, - if^lloo 

- (axo + 6)3 "^^ " =■ ""^^ " ' 

where a < 1. Thus, (|M)) has a unique solution in Sp- Furthermore from ([M)) . it is 
clear that iy(a;) — C as a; — )■ +oo. Recalling the definition of H wc see that the 
asymptotic behavior of g is as given by Lemma [6l | 

Remark 3. From Lemmas [6] and [71 it follows that for given a, h with a > 0, there 
exists a one parameter (C) family of solutions g to (|83p for which g — > as x — > oo. 
Any such solution has the asymptotic behavior given in Lemma [6] 

We seek to prove Proposition[3l For that purpose, for a > 0, recalling the change 
of variable t = t{x) in ([9]), we make the transformation: 

(97) F{x) ^ax + h+ -^9(0 

Note that the change of variable involves the parameters a and 5; c only appears 
in the the solution q{t) as shall be seen shortly. The domain t > T corresponds to 

X > — ^ + \J^^- The change of variables ([97| in ([4]) results in q{t) satisfying 
/no^ 1\ d"^ ( ^ 3 q\dq(\ 3 \ 9^ „ 

Equation ([98)) admits two growing solutions \ft and t corresponding to the freedom 
of changing a and h. The only solution for which — > t — cxd, as noted in 

Lemma[6j corresponds to qit) ^ dr^l'^e^*' for some c. 

From the general theory of representation of solutions by transseries [10] ^ it 
follows that any decaying solution to ([M)) has the following convergent function 
series representation for sufficiently large x: 

(99) g(0 = VrQn(i) , where e '''' 



n — 1 ^ 



where the functions Qn are bounded The equations for Qn are obtained by 



formally plugging in (|99p into (|98p , equating the different powers of ^ and requiring 
that Qn be bounded as < — > oo. Only the equation for Qi is homogeneous while 



(6)a slight modification is needed to accommodate the present ODE which violates a non- 
degeneracy condition on the eigenvalues of the linearization; the changes are minor. Also, 
transseries are used to generate the appropriate ansatz and motivate our choice of go , but play no 
direct role in the proofs. 

('''More precisely, Qn are the Borcl sums of their asymptotic power series. 
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for n > 1, the equation for Qn involves Qj for 1 < j < 71 as a forcing term. 
The associated homogeneous equation does not admit any bounded solution, and 
thus the (5„s are uniquely determined from their equations and the boundedness 
condition. The multiplicative freedom of Qi is equivalent to choice of c and therefore 
without loss of generality, we may assume Qi — > 1 ast— >cx). 

This motivates the choice of the approximation go a-s a truncation of ()99|) (we 
choose to retain two terms in the expansion). To prove that this approximation is 
accurate, we define £ = q — qa and show that £ is small for t > T in an exponentially 
weighted L°° norm. We thus define 



(100) 
where 



£{t)^q{t)-qo{t) 



90 (t) 



ce 



Vt 



■Qi{t) 



(101) 
where 

(102) Qi{t) = 2tIo{t) , where Io{t) := 1 
(103) 

Q^it) = -tIo-tl'^+2tJo , where Jo(t) := l-A/2^e2*erfc(V2i) 



7rte*erfc(\/t) = - 



(l + s)3/2 



ds 



(1 + ^/2)3/2 
as t — )■ 00: furthermore 



ds 



Remark 4. It is clear from (|10ip - ll03p that go(*) 

since, by the change of variables we have g{x) = y ^^°^-^ g(t(.T)), by LemmajBl with 
C ~ \/2a3/2c, we have 

(104) q{t) ^ qo{t) , q'{t) ~ q'o{t) , g"(t) ~ q'^{t) , as t ^ 00 



To analyze the fully nonlinear equation (|98p we write the differential equation 
for £ which follows from (|T00| and (|98l) 



1 



1 

A. 

4^2 



f 

1 3 
^2^2 " 4^3 ' 2t 

where the remainder R = R(t) is given by 
(106) 

qo 

2tJ dt^ 



_3 _go_\ d£^ 

4t2 4i2 J dt 



qojt) 

2^3 



£ = --£" 
2t 



4t2 



R 



d3 



di3 



90 



1 



90 



1 3 

2t ^4i2 " 



\ dqa 






J dt 


I 2^2 


4t^ 



^S,^R4{t) , where ^ 



4^3 



90 



i? 



9i 

4^3 



where 
(107) 



5Qi 
2t 



or 

2t 



13Qi 
4t2 



(108) 



i?4(i) = Q2 



Q2 

"2^ i lt2" 



2t2 + 2^3 
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Using (|106p and (|103|) in (jlOSp we get, after some algebra, 
(109) R3{t)^Jo-tI^-I^ 



(110) i?4(t) = I {II + II - 2/o Jo) + i Jo - \h - ^Jo Jo + \ll 

In the appendix (see equations <\21Q\ and (|225p '). it is shown that < i?3(<) < 
i?3,m < 0.02057 and < Riit) < i?4^™ < 0.0009042 for t > 1.99. Instead of using 
a variation of parameter formula for the third order linear operator on the left of 
(jl05p and turn (jlOSp into an integral equation, we find it convenient to define the 
auxiliary function 

£' £{t) 



(111) hit)^e^[^^t8" ^^.^^3/, 

and analyze the equation for h. 

Remark 5. The choice h is motivated by the observation that 

and thus by (|112p . (|104p and Lemmas [5] and [71 g{x) ~ y^Jt^r) '?(^(^)) ~^ implies 
h{t) as t ^ oo. 



Equation (|99| can be rewritten as 

(113) h' = -'^h + e'B{t)£-^-t'/^e'R, 

where 

nui Rm- '?o(i) , 9^.(i) 'Zo(t) 

The function £{t) can be written in terms of h as follows: 
(115) £{t) = <i/2 /* r s-i/2e-^/i(s)dsdT 



CXD 



We write (|113p in integral form 
(116) 



h{t)^ho{t)-f ^^^^h{T)dT+ f e^B{T)£{T)dT^ f ^^^l^^^^ dT^:J\f[h] (r), 
where 

(117) /io(0 = - /* eW^/^R{T)dT. 

J oo 

We will analyze (|116p to find a unique exponentially decaying ft, (cf. (|118p ). and 
then determine £ from (|115p . 

Remark 6. The functions, qo,£,h and R depend on c as well. For simplicity, our 
notation we will not show this dependence, except when needed. 

We will prove that the operator M defined in (|116p is contractive in some small 
ball in the space H defined as follows: 
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Definition 8. LetH be the Banach space of continuous functions in [T, oo) equipped 
with the weighted norm 



(118) 



\h\\ := supte^*\h{t)\ 



t>T 



We now prove some preliminary results needed in the proof of Proposition |3l 
Lemma 9. Fort in [T, oo) we have 



(119) \\ho\\ < |c|3 -R 



^3 , ?n 



1 



1 



3„ 4|c|e-^ 



2 



3Vr 



where Kz.m o.nd R^^m are upper hounds for \Rz{t)\ and |i?4(t)| m [T, oo). 
Proof. Using (|106p and (|117p we obtain 



(121) \h^{t)\< (|c|V2V'li?3.r 



< 



-R4 



--o,m I 2^13/2 ^^4,)n 

implying (|119p After differentiating (|117p with respect to c, and using (|106p . (|120p 
follows similarly. | 



Lemma 10. 

(122) 



< 



-''\\h\ 



dcS{t) 



< 



-3t| 



9*3/2 



dch\\ 



Proof. We note that 



-l/2„-s 



h{s)ds < -r-3/2||/ij|e-3T 



Using the inequality above it follows that 

^1/2 / ^-1/2 / 5-l/2g-.;j(5)rfg < ^i-3/2g-3t||/^|| 

The bounds on dc£ are obtained in a similar way sinee (|115p implies 
dcE{t;c) = ^1/2 r r-1/2 r s-i/2e-^9,/i(s;c)rfsdr 

I 

Let 

(123) 



Vt 



Q2(T) + 0.0944e-^-i^^^ 



Lemma 11. For t >T, qo (cf. (jlOip ) satisfies the following bounds 

qo 



(124) 
(125) 



< ^ (1 + do\c\) < (1 + do\c\) =: go, 



Vt 



dcqo 



< ^ (1 + 2do|c|) < ^ (1 + 2do|c|) =: qo,c. 
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(126) 



qo 
2t 



< 



Vt 



(127) \0.{,'0-'i} 



< 



Vt 



3|c|e-* 

4^3/2 

, 3|c|e-* 



< 



ce 



T 



3|c|e- 



< 



4T3/2 

3|c|e-^ 
2T3/2 



= : qQ,d.m 



-'■ QOM,c,m 



Proof. Using the integral representation for qo, in the appendix it is shown (see 
(Il99l)-(l200l)) that 



(128) 



(129) 



< 



dcQoit) 



< 



Vt 



1 



2VT 



VT 



{\Q2iT) 



0.09044e" 



-soT 



(|Q2(T)| +0.090446"'*°^) 



where so = 6.159 • • • . implying (|124p and (|125p . Straightforward calculations show 
that 



90 (i) 



-(I-/0) 



1 / 



2t 2t 

where Iq is defined in (|102p and /i — 2tIo. The integral representation (|102p implies 
that Ii := 2tIo e (0, 1) and that 1 - /i = ^/2, where 

hit) = t 

and thus 



2t 

(l + s)5/2 



2t 



V 



ii-h)- 



4t2 



From the equation above and its c-derivative, and the fact that /o,/i,/2 are in 
(0, 1), dull) and ((1271) follow. | 

With do given in (fTM)) let 

(130) d, = Mt-3/2 (1 + |c|do) ; d,,, = (1 + 2|cMo) 



Lemma 12. For i > T > 1, we /ia?;e 



(131) 



g o(T)e^ 
2r 



/l(T)dT 



2t 



h(T)dT 



< Vll^ll : 



Proof. The result follows from Lemma [TTl which implies 



3"go(T) 



2t 

e'^dcqoir) 



< 



2t 



< 



2t3/2 

1 



2t3/2 



(l + |cMo) 



(l + 2|c|fio) 



and noting that \h{t)\ < t-^e~'^*\\h\\. | 
Let 



(132) 



V{t;c) = te*B(t;c) 

c 
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Lemma 13. B defined in pi4|) satisfies the following inequalities for t > T > 1 



(133) 
(134) 
(135) 
(136) 



B{t;c 
dcB{t- c) 
dt 



< 



2t 



< 



2t 



3|c|e~* 

4i3/2 

3|c|e"* 



< 



c e 



< 



e 



2T 

-T 



3|c|e 



-T 



2T 



-B{t;c) 
c 



< e-M 1 



t3/2 



< e 



-T 



4T3/2 

3|c|e-^ 

27-3/2 



— : i?„ 



T3/2 



=: B 



9. 



-B{t;c) 
c 



3e-2* 3e-2^ 

< r- < r- =■■ B. 

- 4^3/2 - 42^3/2 



m,2,c 



V{t; c) defined in (|132p satisfies 
(137) 



(138) 



(139) 



r(t;c) 



Vit;c) 



> 1 



< 



ce 



acF(t; c) 



< 



2T3/2 

4T372 



Proo/. Using (fTOTj) in (fTT4|) it follows that 
(140) 



B{t;c) = - 
(141) 
(142) 
(143) 



2t 8^5/2 (3-^2 - /i) , where /i = 2tla = t 



{l + sf/^ 



ds,h{t)=t 



oa „ — st 



e-'*ds 

(1 + s)3/2 



d,Bit;c)^^ — 



2t 4^5/2 



(3/2-/1) 



2tB{t; c) 



ce 



-2t 



a, 



2tB{t; c) 



4^3/2 

,-2t 



(3/2 +/l) 



r(3/2-/i 



C 4^3/2 ' 

from which, again using the fact that hjli G (0, 1), (|133p - (|136p . To prove (|137|) - 
(1133, note that 



(144) 

(145) 

(146) 
I 



2 cp^* 
V{t; c) --te'B{t- c) = 1 + ^ (3/^ - h) 



V\t-c)^ 



ce 



dcVit- c) 



2t3/2 



/i(0 



4*3/2 



(3/2 - h) 
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Lemma 14. For T > 1 we have 

" <dB\\h\\, 



(147) 

where 
(148) di 



[ e'' B{T)£{T)dT 

J OO 



I e-'dcB{T)£{T)dT 

J OO 



< df 



\h\\ 



(l + |c|rfi), di 



g-T 3g- 

(1 + 2|c|di) , where di 



54T3/2 ' l^l"-!^' "-".c 54^3/2 V" ' ^l^l^-i/' -i 4^3/2 

Proof. Using (jl22p ()133p and (|134p in Lemma [13] the result foUows immediately by 
integration. | 



Lemma 15. For T > we have, 

h{T)£{T) 



(149) 



2r 



dT 



< 



-3T 



90T5/2 



hit) 



-2t\ 



Proof. Using Lemma [TU] and the fact 
fact that h E 7i, the result follows by integration. | 



h\\ (which follows from the 



Proposition 16. For |c| < j, e ~ 0.03 and T > 1.99, there exists a unique 
solution to the integral equation (|116p in a ball of size (1 + £)||/io||i implying that 
ll^ll < (1 + e)||/ioll < 1-6667 x lO"^. 

Proof For T > 1.99, |c| < i and e 0.03, by Lemmas l9l fT2l and [Til and using 
Q2(1.99) = 0.0147- •• and inequalities i?3,„ < 0.0205666, i?4,™ < 0.009042 for 
T > 1.99, see ([216]) and ([225]) in the appendix, we get 

C-3T 



(150) \\Af[h]\\ < \\ha\\ + id, + dB){l + e)\\ho\\ + 



90T5/2 



{i+er\\hor<{i+s)\\h4 



\mhi]-u[h2]\\ <\d,+dB + 



r(l + £)l|/^o|| lift. 



(151) ll-'vL-iJ ■^-r-'-^iU- 1-9 ' I 45^5/2 

implying contractivity of the integral operator in the stated ball. 



,-3T 



(8) 



I 



Remark 7. For |c| < ^ and T > 1.99, heV. implies 
(152) 

\h{t)\ < \\h\\t-'^e-^* < 1.6667 x 10"^ x T-^e"^^ =: /i„ < 1.5651 x 10"*^; Vt > T 

Remark 8. By uniqueness, this is the only solution with /i — > as t — > oo; we 
have proved that such a solution h Cz Ji 

Remark 9. Proposition [TC] extends to any c if T is large enough, as seen in the 
next proposition. This is likely to be useful in extending the present techniques to 
more general initial conditions than ([5]) 

Proposition 17. For any c, there exists T > 1 large enough so that the integral 
equation (jll6p has a unique solution h € Ji. 

Proof. It is clear from Lemmas [9l [121 [T4| that for any given c, the functions dq, ds 
and 1 1 /loll are decreasing in T. Thus, the conditions p50p - (|15ip are met for any 
fixed £ > 0. I 



(8)The values of the error function can be calculated using, for instance, [T], 7.1.28. 
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Lemma 18. For < a < Ur, |c| < ;| and t > T > 1.99, the function £ (see (|115p ) 
satisfies following bounds 



2t 



< ^l.e-''\\h\\ < 1.69 X 10-4-2e-3* 
^ & < ^e-3*t-3/2||/,|| < 9.20 X lO-h-'/'e- 
< V2al/'-\\h\\t-^e-^' < 5.02 x IQ-H-^e 



dx^ 



3t 



-3t 



Proof. From (|115p . the first statement follows immediately. The second statement 
follows from noting that the transformation (|97p implies that 

r-^/'^e-^h{T)dT 



dx V 2i 
Furthermore, we can check 

d2 



da;2 y 2t(a;) 
and hence the third statement follows. | 



Lemma 19. The function h satisfies 



|/i'(-;c)|| <2d,|l/.|i 



Br, 



,-3T 



9T1/2 I 

Proof. We note from (|113[) 



18T5/: 



Hl/^ll 



ce 



1-3, ?7l 



Ti/2 



,-3T 



Using Lemma fT2l and and fT3l we get the statement follows. | 
Lemma 20. The function h satisfies 

-3T 



dch{-,c)\\ < [ 1 - dq - ds - 



457^5/2 



\\h\\ {\\dM\ + idg,c + dB.c)\\h\\oo} 



Proof. We note that pi6p implies 



(153) dch{t;c)^dchoit;c)- J —d,qo{T;c)h{T;c)+l e'dcB{T;c)E{T;c) 



qoir; c)dch{T; c)dT - / e'' B{T;c)dcE{T;c)dT 



2t 

dr 

dch{T;c)E{T;c) — 

Applying Lemmas [H HI [12] and [13] to ([T53]) we get 

(154) 



h{T;c)dcE{T;c) 



dr 
27 



\\dch{t;c)\\ < WdchoW + (dg^c + dB,c)\\h\\ + [dq+dB + 
I 



^-3T 



45T5/2 



\\h\\ \\dM;c)\U 
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Remark 10. Since h d H Lemmas \19\ and \20\ imply 
(155) 



h'{T;c) 

(156) 
dcHT; c 



-2T 



2d,\\h\\ 



M 



-3T 



18^5/2 



R 



3,m 



rpl/2 



■Ra 



1-2T 



1 — dn — du r-77; I 



\\h\\ {\\d,ho\\+{dg,c + dB.c)\\h\\oo} 



3.1. End of proof of Proposition[3l For given |c| < | and T > 1.99, Proposition 
[16] implies that £ satisfies (|99p for t > T > 1.99. This imphes the existence of a 
solution q — qo + £ satisfying (|98p and having i~^/^e~* decay as t — )■ 00. From 
Lemmas [6] and [71 this is the only solution for which as t — > 00. Proposition 

[3| follows since the transformation (|97p for a > in the regime t > T > 1.99 
guarantees that F will satisfy (|4]) for x > |. 



4. Matching and proof of Proposition [4| 

In order for the two representations (|59l) and (|97l) to coincide at a; = | we match 
F and its two derivatives; from (|10ip . (|100p and (|llip we get 



(157) a^F'(^)~a(q'a{Un;c) 
(158) 



go(tm; c) 



— —h(T; c)dT =: Ni(a, b, c) 



b = Fily^N,^^l^qoit^;c)-^l^ 



r-1/2 



e ''h{s\c)ds:=N2{a,b,c) 



(159) 



V2a3/2 



e*^F"{l)^:N3{a,b,c) 



Definition 21. FKe define A = (a, |, f), Ao = (ao, ^, f ) anrfN(A) = (iVi, iiV2, ^iVs) 
Define also 

5a {IIA-A0II2 <Po 5x10-5} 
where ||.||2 is the Euclidean norm and let 



(160) 



^daN2 dbN2 dcN-2 



Note 11. We see that A G Sa implies {a,b,c) G S. The system of equations 
( |-?5?p -| f75P|) is written as 

(161) A = N[A] 

We define J = ^ to be the Jacobian and ||J||2 denotes the P (Euclidean) norm of 
the matrix. We note that 



(162) II J||2 = {daNiY +4.{dbNiy +4.{dcN,y + - {daN2r + {dbN2f + (^c^Va)' 
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Lemma 22. The inequalities 

(163) ||Ao-N[Ao]||2 < (l-a)po 



(164) 



sup ||J||2 < a < 1 
AeSA 



for some a G (0, 1) imply that A = N[A] has a unique solution for A G Sa- 

Proof. The mean- value theorem implies 

(165) 

||N[A]-Ao||2 < ||N[Ao]-Ao||2 + ||N[A]-N[Ao]||2 <Po(l-a) + ||J||2Po <Po 

and also, if Ai, A2 £ Sa- 

||N[Ai] - N[A2]!|2 < IIJII2IIA1 - Aailz < Q'HAi - A2II2 

Thus, (|163p and (|164p imply that N : Sa — > Sa and that it is contractive there; 
the result follows from the contractive mapping theorem. | 

4.1. Proof of Proposition[4l Proposition^] follows from Lemma once we show 
that (|163p and p64p hold. In the following two subsections it will be shown that 
a < 0.764 and that ||Ao - N[Ao]||2 < 1.16 x 10"^ < (1 - a)po thereby completing 
the proof of Proposition 

Remark 12. Note that the proof of Proposition |4] only requires smallness of the 
norms of h and E (we recall that F ^ Fq + E) and on no further details about 
them. If in some application Fq needs to be made C^, then this can be ensured by 
iterating N with h ^ E = 0; the first thirteen digits obtained in this way are given 
in (HH). 

4.2. Bounds on ||N(Ao) - Ao||. We note that 



(166) ao - A^i(ao,6o,co) < ao - ■F'olf ) + «o '?o(im,o; co) - 



qo{tm,0: Co) 



2t 



m,0 



|S'(|)l + yC^V3*"'1 



h\\ < 4.81 X 10" 



(167) 



6o-Af2(ao,&o,co) < fco-i^o(|)+f ao+ 



■9o(^m,o; Co 



Eil) 



9V 2*™'°'^ 



5 

-- ao-iVi(ao, 60, Co) 
< 1.64 X 10"^ 



(168) Co - Af3(ao,&o,co) 



< 



Co 



3/2 



(^(t™,o,co))"'e*"'°F^'(|) 



v2coao 



-2t„ 



\h\\ 



3/2 



^"(f) K(W;co) 



Co^m,0 
.-2t„.o||;,| 



e 



implying 
(169) 



Ao - N(Ao)j|2 < 1.16 X 10"^ 



(^(im,0)"' 



< 1.33 X 10"^ 
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4.3. Bounds on the derivatives of Nj and on | J| 2- We note that 

'25 P 



(170) daNi^-iq'o{t^;c) 



2tr. 



;h{T; c)dT 



-1/2 



2t 



1/2 \ 4 



2 e ''"/i(<„0 



Remark [7] implies that the last two terms on the rhs of (|170p arc bounded by 
C^^e-'-.'/i^ 1 + ^ (f - ^) ■ Applying now Lemmas m [H] to pTUl) we get 



(171) 



daNi 



< 



I'oitmfi, Co) + 



Qoitmfi, Cq) 



1 25 , 

+ i Ty^'r ~ ^o) 



/^25 hi S .I. 



'^'-r ■ 2 



25 

T 



m,r-^rii ^ '^O I . 2/1 1 /9 ' 1/2 



+ qo,d.c.m + a,Q 



25 _ &2 
4 a: 



yWSm.c) (Cr - Co) < 0.081 , 
0/ / 



The maximal value of the bounds is attained when c = Cr and T = tm.i, which we 
used to get the results above. 

(172) dbNi = V2^[2tmBit,n;c) - e-*"/i(t™; c)] 

hence using (|152l) and Lemma [T51 



(173) 



dbNi 



< V2a^ {2tm,rBm + e'^^'^hm) < 0.059 



(174) dcm = -adJq'o{t,n;c)- 



2U 



— a / —j=dch{T; cjdr 



implying from Lemmas 1111 and equation (|156p . 



(175) 

We now consider 
(176) 



dcN, 



< arqo,d,c,7n + are *'"''tJl'^hc,„i < 0.163 



1 



2 V 2ai^ 

1 rr" 



1 r^/25 

2 V 2ir It ^ ^ 



qo{t,n) 

Of 



2 V 2aL 
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It follows from Lemmas O [HI [TH Proposition [TH and equations (|152|) and (|156|) 

that 



(177) 

daN2 



< 



1 



1 



2 V 2ait 



daNi 

25 



1 



1 



2 \/ 2aitm.i 



qo{tm,o', co)+ao 



25 bl 



<Zo(^m,0; Co) 



'ZoCim^o; Co) 



2im.O 



(a^ - ao) + — - 

ai ao 



1 + a„ 



25 bf 



1 

+ 7^ 



1 



2 y 2(iit^ 



90,^,™ + H Oo 



1 



1 



2 y 2aitra,l 



go,. 



ao 



We also note that 



(178) dbN2 = ~-dbNi - (^g[,(t™; c) 

and therefore Lemma [TT] and p52p imply 

5 



25 bl 



qojtm; c) 
If 



qo.d,i 



(cr - Co) < 0.232 



dr 

7^ 



e '^hiT;c), 



(179) 

Now, 
(180) 

dcN2 



dbN, 



< 



dbNi 



-dcqo{tm\c) 



2 V 2t,;_ 

Hence, (|156p and Lemma [TT] imply 



/ s-i/2e-"9c/i(s;c)dsdr 



(181) 
Also, 

(182) daN^ 



dcNo 



< 



1 



dcNi 



2tr. 



-qo,. 



2V2a3/2 

A. 



V^(t,„;c) + -h{tm;c) 



1 -1 



c^'"F"(^) - 



< 0.468 



25 fe2 



l/'(i,„;c) + -/i'(t„,;c) 



Therefore, from Lemma [131 equations ([152]), (|155p and the positivity of Vmin — 

F^' and ^ - ^ - (see Definition[Tl ([26l)-([22l) and ([137]) ): in ([137]) we used c = c^, 
r = tm,i, which minimize Vm)- It follows that 

(183) 



daNs 



< 



2V2( 



3/2 



Cl 



'-(^'<i)-K'i)|){?-|-^(f-S 



25 fe? 



T/ ■ - — 

c; 



Cl 



< 0.44 
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Further, 



(184) 
dbN3 



1 



V{tm]c) + -h{tm;c) 
c 

Using again Lemma [T51 and equations (|152p . (|155p we get 
(185) 



c c 



dbNs 



< 



(f^'(|) + |£;"(|)|) 



14 



} < 0.384 



Furthermore, 
(186) 



V(trn;c) + -h 

c 



^/2a3/2 

Lemma [T3l and equations (|152p . p56p imply 
(187) 



dcV{tm;c) + -dch{Un]c) - \h{tm;c) 



dcN: 



< 



V2( 



3/2 



(F^'(|) + |i?"(|)|) 



^c.m H" hem ~^ 9 

ci cf 



< 0.0029 



By straightforward calculations we get ||Jj|2 < 0.764. 

5. Appendix 

5.1. Bounds on and dcQo- Using the integral representations of /q and Jq, 
PITT)) . (Unil) and ([TOg]) imply 



(188) qo{t) = — e-*t 
where 

(189) U{s) 

Note that 
(190) 



Vt Jo {1 + s)3/2 Jo 



1 



2(1 + 5/2)3/2 2(l + s)3/2 (2 + s)2yTT7 
Q2(t) 



Making the change of variable 
(191) 



s = -1 



4y 



which maps one-to-one (0, oo) onto (3 + 2\/2, oo) we obtain from ()189p 
(192) 
where 

(193) Psiy) = -y^ + (17 + 10^2)^^ - (n + 4V2)y + 3 + 2^2 



TTf ( V/^(3-272)(y-3-272 ^^ 
U^<y)) = (,- 1)3(1 + ^3(.) 
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This cubic has only one real root yo = 30.604- •• implying s = sq = 6.159- ••. 
Since P{y) > for y < yo and P(j/) < for y > yo, so > s > corresponds to 
yo > y > 3 + 2\/2 and s > sq corresponds to y > yo, it follows from (jl92l) that 

(194) U{s) > in (0, so) and U{s) < in(so, oo) 

Further, for s > sq, {i.e. y > yo ~ 30.604..), the functions 

y-^p,iy), ,y-and 



{y-l) {y-i-2V2){y-lY{l + yY 

are decreasing. Therefore, it follows from (|192l) that for s > sq, {y > yo) 
(195) 

(y- 1)-^(1 + y)4 (yg _ 1)3(1 + yg)4 

Therefore, 

0944 

(196) 0> / Uis)e-''ds> ■ e-'°' 

J so * 

Since U{s) is being positive on (0,so) (see (|194l) ) 

(197) f ° U{s)e-'*ds 

Jo 

is clearly a decreasing positive function of t for t > T at thus attains its maximum 
at i = T. Therefore for i > T we have 

poo nSo poa n r\(\AA 

(198) / U{s)e-''dt< / U{s)e-'^ds< / U{s)e-'^ ds + -—e-'°^ 

.In .In .In J 

^2 



We conclude that 



(199) qoit) 



< ^ + ^2^-2^ Q2(T) ^^-2*0 0944g-.oT 

- T T 



We note that that since sq = 6.159- • • , for T > 1.99, 0.0944e-'°^ < 4.5 x lO"'^. 
The numerical value (52(1-99) = 0.0147- - - can be easily obtained from rigorous 
formulas (e.g. [T] 7.1.28) By differentiating qo with respect to c we get in a similar 
way 

P~t ol I -2t -2t 

(200) dMt) < ^ + ^^^Q2(T) + 2|c|^0.0944e-^«^ 

5.2. Bounds on i?3 for t > T > 1.99. In the formula for R3, Iq and Jo 

have integral representations, see ()102p and (|103|) . Using these representations we 
obtain 
(201) 

4 ' io (l + r)3/2(l + (,s_r))3/2 7, (, + 2)2^1 + 7 

(202) 



(.s + 2)2yrT^; 7o 2(s + 2)3(s + 1)3/2 



Adding the expressions above, we obtain from (|109p 

/>oo 

(203) i?3(i) = / dse-^*7^3(s)^^s 
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where 

1 s 3s2 - 4 



(204) 7^3(s) = 
We define 



4(l + s/2)3/2 (.s + 2)2VTT^ 2(s + 2)3(s + 1)3/2 



(205) s = -1 , 

Ay 

Clearly s is increasing in y and maps (3 + 2\/2, oo) to (0, oo). Thus 

(onr\ ( ( \\ 4y3/2(2-V2)(^y-3-2^/2) ^^ ^ 

(206) 7^3(s(y)) = (y + i)6(^_i)3 ^5(y) , 

where 

(207) P^iv) = l-(21-10V2)y+(42-2V2)?;2_(42 + 2V2)?;3 + (10V2 + 21)/-?/^ 

Elementary inequalities imply that P5 has a real root j/o G (33.851,33.852) corre- 
sponding to s = So € (6.9701, 6.9704) By factoring out y — yo and rexpanding the 
rest in powers of y — 5, we find 

(208) P^iy) ^{y- yo) {Ao + Ai{y - 5) + ^2(2/ - 5)^ + A^iy - 5f - {y - 5)^) 
where 

(209) Ao = -y;^ + (10\/2+16)yi] + (48V2+38)?;^ + (238V2+232)?;o + 1139+1200\/2 

(210) Ai = -yl + (11 + IQVTjyl + (98\/2 + 93)yo + 728\/2 + 697 

(211) yl2 = -y^ + (10V2 + 6)yo + 148\/2 + 123 

(212) A3 = 1 + 1072 - yo 

It is readily checked that for y^ in the interval (33.851, 33.852), the coefficients Ao, 
A3 are negative, implying that there has only one zero for ?/ > 3 + 2y/2 > 5, 
namely y = yo- This immediately implies that for s € (0,so) where so = 6.97- • • 
we have 7?.3(s) > 0, while for s > so, Tl-z < 0. 

We now minimize 7^3(s(y)) for y > yo = 33.851 • • • . By simple estimates of the 
derivative, y~^P5{y) is seen to be decreasing, and > y~^P5 > —1 for y G [yo, 00). 
In this interval y — 3 — 2\/2 < y — 1 and thus for s > so, (or y > yo) we obtain 
using ([2061) . that 

V > V'^'(2-V2) Ayl'/\2-V2) 

since 4?/^3/2(-2 _ 2)^^(2/ — 1)^2 jg decreasing for y e (yo,oo). Therefore 

since for t > T, 

/•So /-So 

(214) 0< / e-"*7e3 (s)ds < / e-''^7^3(s)ds , 

Jo "'0 

(215) < - / e-7^3(s) < ~ / e-^' n,{s) < ^I^e'^^^ , 



< - / e-^*7^3(s) < ~ / e-'^TZ^is) < 

Jsa Jso 



T 



it follows that 

(216) R:i[t) < i?3(r) + O-O^O'^ e^soT < ^^(^y-j ^ ^ q2 ^ iQ-s =: i?, „ < 0.02057 
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5.3. Estimating R4. Using (|T02|) -(fT03 l) in (|TTO| we get 

(217) i?4(t) = ^(Ji - h) ^^-h)-^{Ji- ID + f 

where 



(218) h{t) = 2i/o(i) - t ]^ (i + s)^ '^' ^ 



(219) hit) = 4tJo{t) = t ^-—j^ds e (0, 1) 
Clearly Ji{t) > Ii{t). This implies that 

(220) - - ^ ^ ^(^1 - + 1 
Now, 

(l + s)3/2-(l + ,s/2)3/2 



1 /■ 1 - 

(221) V 



(1 + 8/2)3/2(1 + 

which is clearly decreasing in t as is /q (since they are Laplace transforms of positive 
functions), and therefore they attain their maximum at t = T > 1.99. Furthermore, 
it can be checked that 

(222) ^{Ji-If)^^, I e 



2^ _L / „-st 
16*2 V- St Jo 



j 16sVTT^ + lOs^vTTs + 2s^y/TTs + 3sV4 + 2s + 8VTT^ - 8v/l + s/2 ' 
1 (2 + s)='V4T27(l + s)'/' 



Arguing in the same way, we see that ( Ji — /^) is decreasing and attains its 

maximum at t = T, as does Evaluating Iq and Jo (see, once more, [1], 7.1.28) 

we get 

(223) 

< ^(Ji(t)-/i(t))+^ < ^(Ji(r)-/i(r))+:^ < lj„(r)-i/o(T)+:^ < o.oo9042 



and 



1 / r .2W 1 ^ w^x ,2,^^^ _ ^o(T) /2(T) 



(224) < — (Ji - /2) < ^ (Ji(r) - /2(r)) = ^ - < 0.00572 
and therefore, for t > T > 1.99, we have 

(225) Riit) < i?4,m < 0.009042 
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